An attractor for dark matter structures 
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ABSTRACT 

Cosmological simulations of dark matter structures have identified a set of universal profiles, 
and similar characteristics have been seen in non-cosmological simulations. It has therefore been 
speculated whether these profiles of collisionless systems relate to accretion and merger history, 
or if there is an attractor for the dark matter systems. Here we identify such a 1-dimensional 
attractor in the 3-dimensional space spanned by the 2 radial slopes of the density and velocity 
dispersion, and the velocity anisotropy. This attractor effectively removes one degree of freedom 
from the Jeans equation. It also allows us to speculate on a new fluid interpretation for the Jeans 
equation, with an effective polytropic index for the dark matter particles between 1/2 and 3/4. If 
this attractor solution holds for other collisionless structures, then it may hold the key to break 
the mass-anisotropy degeneracy, which presently prevents us from measuring the mass profiles in 
dwarf galaxies uniquely. 

Subject headings: dark matter - galaxies: halos - methods: numerical 



1. Introduction 

Cosmological dark matter structures which 
have been formed in an expanding universe ap- 
pear to have virtually universal profiles. This was 
first emphasized for the d ensity profile through 



is a power-law in radius 



P 



(2) 



numerical simulation when iNavarro et al.l ( 19961) 



demonstrated that all structures from galaxies to 
galaxy clusters are well fitted by the same simple 
form 

P= } ^2. (!) 

J- 1 + X. 



wherc a = 1.875 is in agreement wit h the predic- 
tions from the spherical infall model (jBertschinger 
1985). This numerical evidence was soon used 
in conjunction with the Jeans equation to derive 
density profiles in general agreement with eq. (U) 



dHansenll2004H Austin et al.l2005tlDehnen fe McLaughlin 



where r s is the scale radius. Even though vari- 
ous theoretical attempts have been made to ex - 
plain such density profile (e.g . Henriksen (l2007t) : 
Gonzalez-Casado et all (120071 ) ; lrle fc Kangj (J2010lk 
Hiorth fc Williams! (|2010l )) there is still no con- 
census of the underlying physical reason for this 
universality. 

To ease the task of theoreticians one should 
look for simple connections in some parameter 
space. One such conne ction was suggested by 
Taylor fe Navarro! (|200ll ). who presented numer- 



2005). One problem does, however, present itself, 
namely that recent high resolution simulations 
do not confirm that the pseudo phase-space den- 
sity is a universal straight line ( for recent discus- 
sions, sec Schmid t et al.l (120081 ); iKnollmann et al. 
): lLudiowet all ^OlOh l. 



ical evidence that the pseudo phase-space density 



In this Letter we study the evolution of pure 
dark matter structures using numerical simula- 
tions. We find that all the various equilibrium 
dark matter structures, when perturbed in a sim- 
ple but realistic manner, move towards a particu- 
lar curve in the 3-dimensional space of the deriva- 
tives of density and velocity dispersion, and the 
velocity anisotropy. 



2. Structures in equilibrium 

One of the non-trivial aspects of cosmological 
dark matter structures formed through hierachical 
mergers in an expanding universe is that they are 
seldom in perfect equilibrium: the first moment 
of the collisionless Boltzmann equation, the Jeans 
equation, is not always exactly obeyed unless all 
velocity terms are included. When bulk motion 
can be i gnored , this Jeans equation reads (see e.g. 
iBinnevI (fl982h 



(7- 



2/3), 



(3) 



where we have defined the derivatives of the den- 
sity and radial velocity dispersion, 



dlogp 
dloer 



and ft = 



dlogcr^ 
dloer 



the velocity dispersion anisotropy is given by 



1 



'rad 



and the circular velocity is 

2 GM tot 



(4) 



(5) 



(6) 



Numerical simulations of galaxies or clusters of 
galaxies often have differences between the left and 
right hand sides of eq. Q, of the order 20% for 
large parts of the structures. 

This problem can be avoided by manually set- 
ting up structures in equilibrium. There are sev- 
eral ways of setting up structures which obey 
the Jeans equation. The first is the Eddington 
method, which for a given assumed density pro- 
file creates a unique ergodic structure i.e. de- 
pending only on the energy of t he particles, and 
hence having /3 = by definition (JEddingtonlll916t 
Binnevlll982i ). This method gives the exact distri- 
bution function, /(£), and is therefore often sta- 
ble when evolved by an N-body simulation. When 
we discuss stable configurations, we intend struc- 
tures whose profiles, p(r), of and /3(r), remain un- 
changed when evolved for several dynamical times 
by an N-body code. 

The other method is to make an assumption 
for the radial profiles of the density and veloc- 
ity anisitro py, and then sol ve the Jeans equation 



to get of (JHernquistl 11993). In this method the 



velocity distribution function is assumed to be a 
Gaussian, and the system will therefore not be ex- 
actly stable when evolved by an N-body simula- 
tion. However, after having run for a number of 
dynamical times, such structures settle down in 
an equilibrium configuration, which often has a 
fair resemblance to the original density and ve- 
locity anisotropy profiles. Below we will use both 
methods of defining initial conditions. 

3. Flow towards an attractor 

The Jeans equation in eq. © has infinitely 
many solution: there are 3 unknown functions, 
p(r),af(r) and (3(r), and only one equation. It 
is therefore no mystery that one can create a wide 
range of different structures, with different den- 
sity and velocity profiles, which all obey the Jeans 
equation. The question we are posing here is, if all 
these structures, if perturbed in a simple way, will 
flow to new points in a large volume in solution 
space, or if there is one or several attractors which 
the structure will be drawn towards. 

As solution space it is natural to consider the 
4 dimensional space spanned by r, 7, ft and /3, see 
eq. ([3]). Whereas the Jeans equation allows solu- 
tions to fill a large 3-dimensional volume, then if 
an attractor does exists, we should expect all fi- 
nal structures to land on a low dimensional sheet, 
potentially a line. In practice we divide each struc- 
ture into radial bins. For each radial bin we can 
derive 7, ft and /?, and this radial bin thus rep- 
resents a point in 3d space. In this way we also 
avoid any potential problems concerning scaling 
different structures in radius. 

During the formation of cosmological structures 
in the real universe, these are continuously being 
perturbed through the process of mergers. This 
implies that if an attractor exists in solution space, 
then we might expect repeated perturbations, such 
as the ones from ongoing mergers, to support the 
flow of the structures towards the attractor. 

In order to mimic the perturbations of the 
mergers in a controlled way, we perturb the equi- 
librated structure in a simple yet fairly realistic 
manner. The velocity of each particle is multiplied 
by a random number, in such a way that any given 
radial bin has exact conservation of kinetic energy 
and angular momentum. We assure that no parti- 
cles get a velocity above the escape velocity. Natu- 



rally, the potential energy is also conserved during 
the perturbations. Now, having perturbed each 
particle in the structure in this controlled man- 
ner, we start the N-body code, and let the struc- 
ture settle into a new equilibrium condition. We 
let each simulation run for a time corresponding 
to one dynamical time, tdyn = ^1 \fG~P~, at radius 
13 times the scale radius. 

The simulations were performed with the public 
versio n of Gadget-2 (|Springelll2005t ISpringel et al. 



20011 ) which is a massively parallel N-body code 
that uses a hierarchical tree algorithm to calculate 
gravitational forces, and gives individual timesteps 
to all particles. A spline softening of 0.005r s was 
used. We are being overly conservative in the 
analysis, and include only regions outside of 5-10 
times this softening. We also require that the bins 
contain at least 5000 particles. 

4. The wide range of initial conditions 

In order to allow the systems to start from 
a large range of initial conditions, we initiate a 
set of significantly different systems. Using the 
Eddington method, we have an isotropic system, 
whic h has a density profile like a Hernquist pro- 
file (|Hernquist]ll990l ). with inner and outer density 
slopes of — 1 and —4. We construct a cored profile, 
with central slope of 0, and outer slope of —5. We 
also use a hoovering profile, p ~ l/(l + (r/ri) 2 (1 + 
r/r2) 3 ), where r\ — 0.01 and r^ = 100 which is 
characterized by having a large part of the cen- 
tral structure where the slope is —2. We also cre- 
ate a set of initial conditions with Osipk ov-Merrit 
anisotropies (|Osipkovll979HMerritdll985h with the 
property /3(r) = 1/(1 + (r a /r) 2 ): a Hernquist den- 
sity profile (slope —1,-4) with r a = 1, and a 
cored profile (slope 0,-3.5) with r a — 1.5. Using 
the gaussian method, we create 4 structures, with 
both trivial and highly non-trivial j3{r) profiles, 
all with a Hernquist density profile (slope —1, —4). 
All structures presented here contain 10 6 DM par- 
ticles, and low resolution runs with 10 5 particles 
establish that this is sufficient for our purpose. All 
of these initial conditions (after having run for 3 
dynamical time at radius 13 times the scale ra- 
dius) are presented in figure [TJ This figure shows 
that our initial conditions span a large part of the 
solution space of (/3, 7 + k). The Jeans equation 
allows all of this parameter space to be populated, 



however, there are no initial conditions in the up- 
per right corner, w here the inequality f3 < —7/2 
(JAn fc Ev ans 200 6]) prevent s any stable solutions 
fsee also lCiotti fc MorgantH (J2010I) ). 
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Fig. 1. — The initial conditions. A 2-dimensional 
projection in solution space to the Jeans equation. 
According to the Jeans equation, eq. ([3]) , all points 
in this space may be populated. No stable solu- 
tions exist in the upper right corner. Our various 
initial conditions of equilibrated structure are cre- 
ated to span a large volume in the 3-dim solution 
space. All the structures of this figure are in equi- 
librium, and are stable when evolved in an N-body 
simulation. 



5. The perturbations 

All perturbations follow the repeated pattern of 
kick-flow. First the structure is perturbed (kicked) 
and then the N-body simulation allows it to settle 
into (flow towards) a new equlibrium configura- 
tion. This kick-flow is repeated numerous times. 

We first consider 2 structures, each with ini- 
tial Hernquist profile and with ,9 = 0. These 
are created using the Eddington and the Gaussian 
method respective. We then perturb the energies 
of each individual particle in the range [0.75, 1.25], 
in such a way that there is energy conservation at 
each radial bin. The structure is then evolved with 
Gadget-2 for 1 to 3 dynamical time at 13 times 
the scale radius. This process is repeated approx- 
imately 20 times. The final structures are com- 
pared, and we find no significant difference when 
plotted in the plane (/3, 7 + «). This shows that 



there is virtually no dependence in the final result 
on the method of creating the initial conditions 
(Eddington vs. Gaussian). 

We then take the isotropic Hernquist struc- 
ture, using the Eddington method, and perturb 
by changing energies in the much larger range 
[0.25, 1.75]. Also this structure ends in the same 
region in (/3, 7 + re) after repeated perturbations. 
This demonstrates that the magnitude of the per- 
turbation is not crucial to the final structure. We 
also test perturbations where we don't prevent ve- 
locities above the escape velocity: for these struc- 
tures some particles do leave the structure, and 
these perturbations are therefore less conservative. 
Also structures perturbed in this way land on the 
same curve in solution space. We can therefore 
perturb all structures by changing energies in the 
range [0.25,1.75], irrespective of the method of 
creating the initial conditions. 

The result of all these perturbations is pre- 
sented in figure [2] from which it is clear, that all 
the structures end up along a 1-dimensional curve 
in this 2-dim parameter space. It is important 
to emphasize that we had initial conditions both 
above and below this attractor solution. When we 
plot in the 3-dim space spanned by 7, re, j3 we still 
find that the attractor solution spans a non-trivial 
1-dim curve. 

This curve is clearly S-shaped, but we can ap- 
proximate it with a piecewise straight line, which 
for small (3 could be (7 + re) = —8/3, and for inter- 
mediate j3 it looks like (7+ re) = —0.7—4/3. It is in- 
teresting to note t hat the anisotropy- s lope c onnec- 
tion suggested by lllansen fc Moord (|2006l ) looks 
like (7) = —0.8 — 5/3 in the reg ion which can be 
trusted, namely —2.5 < t < — 1 (IHansen fc Stadel 



20061: IZait et al.ll2008t [Navarro et al.ll2010h . which 
is in good agreement with our findings. For very 
high (3 we are in the very outer region of the struc- 
tures, where full equilibrium can be questioned, 
since a larger spread in the final points is observed. 
The included points have been allowed to equli- 
brate at least 20 dynamical times in total, how- 
ever, the outermost points have only been allowed 
1 dynamical times since the last perturbation. We 
therefore won't speculate on the form of the at- 
tractor there. 

Whereas the 1-dimensional attractor has a non- 
trivial S-shape in the space of 7, re and /3, then 
it happens to lie in a simply parametrized sheet, 
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Fig. 2. — The attractor solution. The same 2- 
dimensional projection in solution space of the 
Jeans equation, eq. (|3]), as presented in Fig. 1. 
The points represent the final profiles for all the 
initial conditions, after these have been perturbed 
and allowed to equilibrate repeatedly. This attrac- 
tor solution is also a 1-dim curve in the full 3-dim 
space spanned by 7, re and (3. 



given by one single free parameter, S — 0.15, 



Y 



(i + Y 3 y 



(7) 



where Y = —Sj — (1 — S)k. This means that when 
seen in a slightly different projection, the attractor 
in fig. ([2]) looks like a very simple curve. 

With this attractor we can effectively remove 
one degree of freedom from the Jeans equation: 
for a given density profile there is now a unique 
connection between re and /3, which allows us to 
calculate both. 

6. An effective equation of state 

The Jeans equation describing collisionless par- 
ticles is fundamentally different from the Eulcr 
equation for collisional fluids, however, a stan- 
dard mapping between the two is made by setting 
(3 = and interpreting m all 2 as a temperature 
( Binnev fc Tremaindll987 ). The non-trivial part 
of this interpretations is to assume that (3 = 0, 
since numerical simulations for yea rs have demon 



strated that (3 is not generally zero (jCarlberg et al 
1997tlColefcLacevlll996l) . 



We now allow ourselves to use our newly found 
connection to make a different interpretation of 
the Jeans equation from a fluid perspective. 

Let us first consider the region with small /3, 
where the attractor solution approximately reads 
(7 + ft) = —8/3. When we insert this into the Jeans 
equation, we get 



GM U 



-a 2 A (7 + n) , 



(8) 



where A = 3/4. If we attempt to describe the dark 
matter by an effective equation of state (general- 
izing a polytropic index) we have 



A'r, 



P = b(pal)"\ 



(9) 



where b is a constant. Thus, if we require that the 
Jeans equation, eq. ©, should take the form of 
the equation of hydrostatic equilibrium, GM/r 2 = 
VP/p, then we find the solution K = 3/4. We 
can therefore interpret this by the DM behaving 
like a fluid with an effective polytropic index of 
3/4. 

For intermediate /3 we found approximately the 
connection (7 + k) — —0.7 — 4/3. This is almost 
solved in the same way, finding an effective poly- 
tropic index of 1/2, however, the constant term 
(B = —0.7) changes the picture. In order to solve 
for such an effective polytropic index, we gener- 
alize eq. © to P = b p Kl a 2K ' 1 , where & is a new 
constant. We now need to introduce an assump- 
tion about the pseudo phase-space density, namely 
pja% ~ r~ a . Doing this we find the solution, 
K x = 1/2 - B/a, and K 2 = 1/2 + Be/2a. 

7. Will any perturbation do? 

A ball will always run downhill, if kicked ran- 
domly by a group of children playing, however, if 
they consistently kick the ball uphill, then the ball 
may appear to defy the attractiveness of gravity. 
We should expect the same to be possible with 
the abovefound attractor. In order to quantify 
this idea we tried to perturb only the radial ve- 
locity component. When the radial component is 
perturbed by a random number (always obeying 
exact energy conservation in any radial bin) then 
the radial velocity distribution gets flattcnd out. 
When we subsequently let the system evolve by 
the N-body simulation, the high energy tail of the 
distribution essentially gets cut off, and thereby 



the effective radial dispersion gets slightly smaller. 
When we repeatedly perform such kick-flow, we 
even manage to get negative values of /3 at inter- 
mediate radii. This confirms our expectation that 
perturbations may be invented, which will defy the 
attractor. If we instead only perturb the tangen- 
tial velocities, then the system flows towards the 
attractor. 

We see from this discussion that it is relevant to 
consider which kind of perturbations should hap- 
pen during structure formation. During merging 
there will be violent relaxation. Two salient fea- 
tures of violent relaxation is that it changes the en- 
ergy of individual parti cles on a timescale similar 
to the dynamical time (|Lvnden-Belllll967 ). which 
for mergers is of the order R/v, where R is the size 
of the smaller merging object and v is its velocity. 
We are here interested in the energy exchange with 
specific particles in the larger structure, which has 
a different timescale from the equilibration of the 
smaller structure. The other feature is that violent 
relaxation acts on all particles irrespective of their 
orientation. In our example above we have an in- 
stantaneous exchange of energies between the par- 
ticles (reminicent of violent relaxation) and then 
we let the system flow (which provides the needed 
phase mixing). 

In our particular setup, the kicks do not change 
any properties of the profiles by themselves, and 
without the kicks there is no effect of the flow. 
Thus, when interpreted as violent relaxation (kick) 
and phase mixing (flow), we see that both are 
needed to move towards the attractor. 

Antonov's laws of stability state that e.g. an 
isotropic Hernquist structure will remain isotropic 
when exposed to minor perturbations. These laws 
are derived under the assumption that the r.h.s. 
of the Boltzmann equation is zero. Our perturba- 
tions act like an instantaneous non-zero term on 
the r.h.s. of the Boltzmann equation, and we are 
therefore not violating Antonov's laws of stability. 

Finally, it is worth emphasizing that we find 
the attractor in the space of 7, K and /3. Very 
specifically, we do not in this work identify one 
unique density profile, p(r). This may be be- 
cause we are constraining the perturbations to al- 
low very little freedom in modifying the density 
profile, or it ma y be because th e universa lity of the 
densi ty profile ([Navarro et al.l 4996; M oore et al.l 
Il998l ) really finds its origin in accretion history 



( Gonzalez-Casado et al.ll2007 ) 



8. Breaking the mass velocity anisotropy 
degeneracy 

Let us finally discuss what our newly found at- 
tractor may do for the mass-velocity anisotropy 
degeneracy. When we observe the stellar kinema- 
tics in a dwarf galaxy we can observe the stellar 
density and the stellar dispersion. Then, the jeans 
equation, eq. ©, tells us that for any assumed ve- 
locity anisotropy profile for the stars, /J(r), we can 
solve for the total gravitating mass. However, if 
we had assumed a different /9(r) then w e would 
have found a different total mass profile (jStrigari 
2010h . 

If we for instance consider a Hernquist density 
profile with a /3 profile in agree ment with numeri- 
cal simulations and obse r vations (Hansen & Piffa rettil 
20071 : iHost et all 120091: IWoitak fc Lokasl l2010h . " 
then the reconstructed mass is overestimated by 
up to 40%, if the analysis is made under the sim- 
plifying assumption /3 = 0. Also the derived inner 
density slope (from the total mass) is systemati- 
cally found to be more shallow than the true slope 
is, by up to 10%. This means that if the true 
density slope is —1, then we will measure around 
—0.95, if we assumed /3 = in the analysis. 

This may no longer have to be the case. If 
our attractor solutions also applies to stellar sys- 
tems in a dwarf galaxy, or to the dynamics of the 
galaxies in a galaxy cluster, then we have a unique 
connection between the 3 quantities, 7, k and /3. 
Therefore, if we have measured (accurately) the 
stellar density and dispersion profiles, then we do 
in principle know exactly what /3(r) looks like, and 
we can then deduce the unique total mass profile. 

9. Conclusions 

We have identified an attractor solution for 
dark matter structures. This implies that any 
dark matter structure which is repeatedly per- 
turbed (e.g. through violent relaxtion during 
merging) and then allowed to relax (phase mix), 
will flow towards this 1-dimensional curve in the 3- 
dimensional space spanned by the 2 radial deriva- 
tives of the density and velocity dispersion, and 
the velocity anisotropy. 

This finding provides strong support for the 



idea that the universalities found in cosmological 
dark matter structures are a property of gravity, 
and not simply a result of similar accretion and 
merger histories of different structures. 

This attractor solution effectively removes one 
degree of freedom from the Jeans equation, giving 
hope that we will eventually be able to solve the 
Jeans equations analytically, and thereby truely 
understand the origin of the universal profiles. 
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